Dissipative hydrodynamics for relativistic multi-component systems 
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Second-order dissipative hydrodynamic equations for each component of a multi-component sys- 
tem are derived using the entropy principle. The shear viscosity of the whole system, appearing in 
the equation summed-up over all components, is related to the partial shear pressures and cannot 
be considered as an external parameter. We demonstrate that it is essential to solve hydrodynamic 
equations for each component, instead of treating a mixture as an effective one-component system 
with a free parameter rj/s. Thus, extractions of the rj/s value of the QGP at RHIC and LHC have 
to be reexamined. 
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The deconfined state of QCD matter produced at 
the early stage of ultrarelativistic heavy-ion collisions at 
RHIC and LHC is a multi-component system with quark 
and gluon degrees of freedom. Large values of elliptic 
flow coefficient V2, observed in experiments of ultrarela- 
tivistic heavy- ion collisions at RHIC P, Q and LHC Q , 
indicate that the produced quark-gluon plasma (QGP) 
is a nearly perfect fluid. This has motivated rapid de- 
velopments on relativistic dissipative hydrodynamic for- 
malisms [iHsl . The value of the shear viscosity to the 
entropy density ratio rj/s for the QGP at RHIC and 
LHC was extracted from comparisons of hydrodynamic 
J| as well as kinetic transport calculations with ex- 
perimental data. All these hydrodynamic formalisms are 
based on the assumption that the quark-gluon mixture 
can be regarded as an effective one-component system, 
where rj/s is an external parameter characterizing the 
dissipation in the system. One may ask the question 
whether this assumption really holds in general, which is 
of major interest for investigation of the QGP properties. 

In this Letter we demonstrate that a standard one- 
component hydrodynamic description with a single shear 
viscosity coefficient in general cannot be applied to a 
multi-component system. We will explain this statement 
by deriving second-order dissipative hydrodynamic equa- 
tions for a multi-component system from the entropy 
principle. Our approach differs from the one reported in 
Ref. [8( , since we introduce separate evolution equations 
and transport coefficients for each component of the mix- 
ture. We then show that by summing-up equations for 
all components one can obtain an equation for the system 
as a whole, which has a relaxation- type form characteris- 
tic for all the present hydrodynamic formalisms, but the 
effective shear viscosity for the mixture is now related to 
the partial shear pressures of its components and thus 
cannot be considered as an externally specified parame- 
ter. It is essential to solve the hydrodynamic equations 
for each component. We will confirm our findings by 
comparing the solutions of the hydrodynamic equations 



with those from kinetic transport calculations. 

We consider a mixture of N particle species. Neglect- 
ing bulk pressure and heat flow we construct the total 
entropy current as [ill ] 
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where s eq is the total entropy density in local equilib- 
rium and m m is the hydrodynamic velocity. 7j and are 
the temperature and local energy density of the particle 
species i. ft = (9/4^) and Trf = T? v - T^ q is the 
shear stress tensor, which is the difference between the 
energy-momentum tensor Tf v and the equilibrium one. 
Equation (JXJ) is the generalization of the entropy current 
for a one-component system (JV = 1). 
The total entropy production is then 
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and A a p = g a p — u a up with the metric g a p = 
diag(l, — 1, — 1, — 1). We have used conservation of the 
partial particle flows and total energy-momentum ten- 
sor, <9 M 7Vf = and d^T^ = 0, to obtain Eq. ©. 

According to the second law of thermodynamics the 
entropy production is non-negative. A simple way to ful- 
fill this is to make the terms in the square bracket in Eq. 



© to be proportional to 7rf p , [•••] = < p / '(2^7*) . Then 
the entropy production © has the following algebraic 
structure: 
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This leads to the dynamical equation for each 7r™^ 
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which was introduced by Israel and Stewart for a one- 
component system (N — 1) 0, 0j- The main subject of 
this Letter is to confirm the validity of Eq. (j4]) for a 
multi-component system (N > 1). Before we do so, we 
shall first determine the coefficients rji, which in general 
differ from the usual definition of the shear viscosity. The 
reason for this is that 7r°^'s are correlated due to inter- 
actions between particles from different species. These 
correlations between tt^'s can only be seen, when each 
r\i depends on all tt"' 3 , j = 1, 2, • • • , N. We will also show 
later that the rji become the shear viscosities, when the 
ratios of components of tt? s are relaxing to constants 
in time. 

We now make use of relativistic kinetic theory and ex- 
press the entropy current via the particle distribution 
function in phase space fi(x,pi): 

N 

s " = E / dT * rfM x >Pi) f 1 - k AC*, ft)] (5) 

with dTj = d 3 pi/ Ei/ (2ir) 3 . It was shown for the case of 
N = 1 [lJ and is obviously true for TV > 1 that using the 
Grad's ansatz ^ fi(x,pi) = fi,eq(jx,Pi){l + Aiir itliI ,pPp?) 
in Eq. ([5]) one obtains Eq. (fTJ up to second order in 7rj )M „. 
Here fi,eq{%,Pi) is the equilibrium distribution function 
and Ai = [2(e^ + Pj)T?] _1 , where Pj is the pressure. 

Assuming that space-time evolutions of fi(x,pi) obey 
the Boltzmann equations 
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where Cjj are the collision terms describing interactions 
of particles of same species and CV, describing binary 
interactions of particles of different species. Explicit ex- 
pressions for the collision terms can be found for example 
in [l3| | . Taking derivative of ([5]) and using (j6|) we obtain 
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Comparison between Eqs. Q and Q leads to 
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Because the collision term Ci is a functional of all //s, 
each rji depends on all ^"'s with j = 1, 2, • • • , iV. 

Equations ((4]) together with ([8]) are the main findings 
in this work. They provide the dynamic evolution of the 
shear pressure for each component of the mixture. This 



is highly relevant for the present research on high energy 
nucleus- nucleus collisions at RHIC and LHC, where the 
produced quark-gluon plasma is a mixture, and thermal 
equilibration of quarks and gluons is expected to proceed 
differently during the expansion. For a one-component 
system, the shear viscosity 77, or practically the shear vis- 
cosity to the entropy density ratio 77/ s, can be regarded 
as a to be specified but free external parameter. In con- 
trast to that, for a multi-component system we recognize 
that each transport coefficient rji depends on shear pres- 
sures of all components and, thus, is not a free parameter. 
This is a new property, which is not observed for a one- 
component system, although the basis of the derivation of 
Eqs. (HJ) is the same as in the second-order Israel-Stewart 
theory for a one-component system, namely the entropy 
principle Eq. ([3]). One may ask the question whether it 
is still correct to regard a multi-component mixture as an 
effective one-component system with a freely chosen 77/s, 
as mostly done in Refs. [9j for the quark-gluon plasma. 

In the rest of this Letter we will discuss some applica- 
tions of the obtained equations. To confirm our findings, 
we solve the multi-component hydrodynamic equations 
and compare the results with those from kinetic trans- 
port theory. 

For reasons of simplification we consider a one- 
dimensional expansion with the boost-invariance. In this 
case w/ 1 = (t, 0, 0, z)/t, where r = yjt 2 — z 2 is the proper 
time. Because of the transversal isotropy the shear stress 
tensor is then 7H 1 = diag(Q,%i/2,'Ki/2,—'Ki) in the co- 
moving frame, iti denotes the shear pressure of species i. 
The hydrodynamic equations (j4]) are reduced to 
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TTj is the derivative with respect to r. The set of equa- 
tions is the multi-component analogof the known Israel- 
Stewart second-order equations 12 , 14 1 and has a 
relaxation-type form, existence of which was advocated 
in Ref. [H)]. 

Furthermore, we consider only binary elastic collisions, 
which keep the particle number of each species conserved. 
Particles are assumed to be massless Boltzmann parti- 
cles. We also take an isotropic distribution for the colli- 
sion angle. With these simplifications we can analytically 
calculate the integrals in Eq. (JSJ and obtain (for details 
see [H) 
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where n^s are the local particle number densities, ct^ 
denotes the cross section for a collision between a particle 
of species i and a particle of species j, and an is the 
cross section for a collision of two identical particles of 
species i. One can clearly see that rji explicitly depends 
not only on temperature, cross sections, and the chemical 



3 



composition, but also on the ratios of the shear pressures 
7Vj / 7r<i . Even when tt^s become small compared with e^'s, 
the ratios TTj/wi cannot vanish, rji is a coefficient rather 
than shear viscosity of the species i in the usual sense. 
We now want to address the question whether and how 
dissipation in a multi-component fluid as a whole can be 
described by a single effective shear viscosity coefficient. 

Summing up the equations Q over i, we obtain the 
hydrodynamic equation for the total shear pressure 7r = 
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with the definition of the effective shear viscosity for a 
mixture via 
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For the second identity the relation Tj = ej/(3rii) and 
Eq. (ITU)) were used, e is the total local energy density. 
We recognize the equality of the two sums in (IT21 and 
thus can write 



Veff = 



where X m fp,i = 1/ Si=i "i "^' is the mean free path of 
particles of species z and = 7Ti /tt is the fraction of the 
total shear pressure. For equal mean free paths, \ m fp.i = 
A, a multi-component system as a whole behaves as if all 
particles were identical. For this special case one finds 
r] e ff — (2/5)(e/A), which is, indeed, exactly the result 
for a one-component system [3, EH, as it should be. 
For unequal mean free paths the effective shear viscosity 
additionally depends on the dynamic values a*. This 
means that in general, it is not possible to describe a 
mixture by a free parameter r) e ff/s. One has to solve the 
equations ((3]) for each shear pressure, rj/s is no longer a 
characteristic value. Mean free paths and temperatures 
of all particle species are the relevant scales. 

On the other hand, the shear viscosity r\ of a mixture 
can be calculated by means of the Green-Kubo formula 
17| and no dependence of r\ on the ratios a^s is expected. 
However, the Green-Kubo formula is applicable when the 
particle system is in thermal equilibrium. Equation (|13[) 
gives a more general expression for the case the system 
is out of equilibrium. When the system is approaching 
equilibrium, a,'s become constants, and their values can 
be obtained by solving dj = using Eqs. ©. The calcu- 
lation is difficult, because we are dealing with a system 
of non-linear equations. We can solve it for the simple 
case of N = 2 and obtain: 
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Putting (fT4")l into (TT3")) gives the shear viscosity of a two- 
component mixture in equilibrium, rj eq . We would like 
to mention that the weights in Eq. (1131) depend not 
only on the composition ni/712, but also on the difference 
between the inverses of the mean free paths, which is not 
trivial. 

To verify our findings we now compare the solu- 
tions of hydrodynamic equations with those from kinetic 
transport calculations using the partonic cascade model 
BAMPS [l^| . Previously published works demonstrated 
that BAMPS results can be reg arded as a benchmark for 
hydrodynamic calculations 12. 18|. 

As a further simplification we consider a spatially ho- 
mogeneous two-component particle system. This means 
that terms containing the gradient d^u^ 1 = 1/r in Eqs. 
© and (fTTj) are dropped. Numerically this is realized by 
confining particles in a static box. 

Initial condition is chosen as follows: T\ = T2 = 400 
MeV, ni/712 = 5, TTi/ei = 7r 2 /e 2 = 0.3. Particles of 
species 1 are in chemical equilibrium, whereas particles 
of species 2 are undersaturated. Cross sections are G\\ = 
3.88 mb, 022 = en/4, and a\% = crn/2. We thus mimic 
elastic interactions among gluons (species 1) and quarks 
(species 2). 

Figure [TJa) shows the time evolutions of the partial 
shear pressures tt\ and 7T2, obtained by solving the hydro- 
dynamic equations © (lines) and from BAMPS (sym- 
bols). Results are divided by total energy density. We see 
perfect agreement between the two approaches over up 
to 4 orders of magnitude. This demonstrates the valid- 
ity of the new multi-component hydrodynamic equations 
derived here. To achieve such high numerical accuracy 
in BAMPS calculations we set the box to have a size of 
16 x 16 x 16 fm 3 and use a test particle number of 2 • 10 2 . 
The total particle number in one run is 1.3 -10 7 , and the 
results shown are the average obtained from 10 4 runs. 

Results for the total shear pressure divided by total 
energy density are also shown in Fig. [JJa), compared to 
solution of Eq. (fTTj) with the shear viscosity in thermal 
equilibrium rj eq . We see that the effective one-component 
hydrodynamic description underestimates the thermal 
equilibration at early times, while it can be applied at 
late times when the system is approaching equilibrium. 
In order to meet the BAMPS results at both initial and 
final times, one has to choose a smaller r\ eq . This gives 
rise to conclusion that the extraction of the rj/s value for 
the quark-gluon plasma at RHIC should be reexamined. 

Finally, the m / 7r 2 ratio is depicted in Fig. [JJb) . Start- 
ing from the initial ■k\/-K2 = 5, the ratio indeed relaxes 
to the value calculated analytically, 7Ti/7r 2 = 0.488. [see 
Eq. I[14p] . The effective relaxation time for tti /7r 2 in this 
case is approximately 0.6 fm/c and has a complicated 
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FIG. 1. (Color online) (a) Time evolution of the partial and 
total shear pressures normalized by total energy density and 
(b) ratio of the partial shear pressures from BAMPS (sym- 
bols) and hydrodynamic calculations (lines). Solid grey lines 
represent the effective one-component solution (see text). 

dependence on densities and cross sections. 

In this Letter we have derived novel hydrodynamic 
equations for each component of a multi-component 
system. Although the equations possess the relaxation 
form similar to that for a one-component system, 
they do reveal a new property - the shear viscosity 
coefficients for each component as well as the effective 
shear viscosity of the mixture depend on the ratios of 
shear pressures. This leads to the conclusion that hydro- 
dynamic behaviour of a multi-component system as a 
whole in general cannot be described by one-component 
hydrodynamic equations with an effective shear viscosity 
coefficient calculated using standard methods (e.g. 
Green-Kubo relations). It is inevitable to solve the 
hydrodynamic equations for each component. Instead of 
T]/s, temperatures and mean free paths are the relevant 
scales. We have confirmed our findings by comparing the 
solutions of the derived hydrodynamic equations with 
those from kinetic transport calculations using BAMPS. 
Both results agree with each other with high accuracy. 
Although the confirmation has been demonstrated for 
a static medium with simplifications on cross sections, 
our findings are expected to hold for general cases. A 
great potential of the present work lies in investigations 
of hydrodynamic behaviour of the quark gluon plasma 
at RHIC and LHC 0. 
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